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Abstract. Various phase shifters and absorbers can be put into the arms of a double 
loop neutron interferometer. The mean intensity levels of the forward and diffracted 
beams behind an empty four plate interferometer of this type have been calculated. It 
is shown that the intensities in the forward and diffracted direction can be made equal 
using certain absorbers. In this case the interferometer can be regarded as a 50/50 
beam splitter. Furthermore the visibilities of single and double loop interferometers 
are compared to each other by varying the transmission in the first loop using different 
absorbers. It can be shown that the visibility becomes exactly 1 using a phase shifter 
in the second loop. In this case the phase shifter in the second loop must be strongly 
correlated to the transmission coefficient of the absorber in the first loop. Using such 
a device homodyne-like measurements of very weak signals should become possible. 
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1. Introduction 

Neutron interferometry is already a well known method for investigation of coherent 
matter wave properties. Many features of neutron waves have been analyzed such as co- 
herence and post - selection effects, the spinor symmetry, the magnetic Josephson effect, 
the multi - photon exchange interaction, spin superposition problems, the Aharonov 
- Bohm effect, topological phases, gravitationally induced quantum interference, the 
Sagnac effect, the neutron Fizeau effect et cetera [1]. 

In this paper we concentrate on coherence effects in a double loop interferometer [2], [3]. 
Fig.(l) describes an extension of a standard three - plate neutron interferometer where 
a second loop is added by using a second mirror crystal M2 and where various phase 
shifters A and absorbers a can be inserted in each loop. Both loops are coupled via 
beam (d) and therefore most attention will be given to the action of an absorbing phase 
shift in this beam. We focus on two phase shifters and A/ and two absorbers ad and 
af in beam (of) and (/) respectively because phase shifters or absorbers in beam (6) are 
of no additional significance. 

The paper is organized as follows. In section 2 intensity aspects of the double loop 
interferometer are considered. In spite of the non symmetric diffraction properties of 
neutrons in single crystals it can be shown that a 50/50 intensity split-up can be reached 
using an absorber in loop B. One of the main findings of the analysis of the double 
loop interferometer in section 3 is that the visibility of K in the forward beam can be 
made very high (ideally 1), even for very weak input signals. In section 4 the results are 
discussed in detail taking into account real experimental situations. 

The main aim is to find an interaction where small signals transmitted through phase 
shifter A^ can be detected with high precision. In this respect it is a search for 
homodyne-like detection of weak neutron signals by the constrain that a symmetric 
beam splitter does not exist in case of diffraction from a crystal. The new system con- 
sists of a coupled double loop perfect crystal system and provides a symmetric beam 
splitting or can serve as a basic of homodyne neutron detection in a similar sense as 
known for photon beams [4]. 

2. Intensity aspects 

From the dynamical theory of diffraction the wave function behind the analyzer crystal 
in the forward direction (index 0) which is a superposition of three waves due to the 
three pathways (s. Fig.(l)) reads as ([1],[5]) 

^acf)+(ad 9 )+(be 9 ) = ^vy^^kAf-a, + e ikA d ~a d + ^ _ ^ 

In this expression ip e = -u elkr is the incoming plane wave function with amplitude Uq, 
wave vector k and position vector r impinging on the beam splitter S. e n is a phase 
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Figure 1. Double loop neutron interferometer: there are three beam pathes: (acf), 
(adg) and (beg). Ad and A/ are two phase shifters, ay and a/ are two absorbers, K 
and Kg denote the intensity in the forward and diffracted direction respectively. S is 
the beam splitter crystal, Ml and M2 are two mirror crystals and A is the analyzer 
crystal. All four crystal plates have the same thickness D and have equal distance 
from each other. The first loop (loop A) is formed by the beams (adeb), the second 
loop (loop B) by the beams (cfgd). Beam (g) is a superposition of beam (ad) and 
beam (be). 



factor which is of no relevance in the following. The crystal functions Vo, vq and v' are 
given as [5] 

.„„ = .""wvn^) + •« sm{A / y^ ) ] . (2) 

v i + y 



. iPD sin{Ay/l + y 2 ) , 
v g = -ie /r-i— 7 > v o = v G (-y) • (3) 

The index G means the diffracted direction. P = — 7r(l +y)/A and A = ttD/Aq where 
A is the Pendelloesung length which is A = A cos{ r y)/\V{G)/ E\. A = 27r/A; is the wave 
length and 7 is the angle between the beam direction and the vector perpendicular to 
the crystal surface. |l / (G)/i?|~10" 6 is the ratio between the crystal potential V(G) 
and the energy E of the particles where G is the reciprocal lattice vector corresponding 
to the reflection plane under consideration. A = 0.00641cm for the (2,2,0)— reflec- 
tion of 2A neutrons in a silicon single crystal. The quantity A is proportional to the 
crystal thickness D. The quantity y depends on the deviation from the Bragg angle. 
For symmetric diffraction y is defined as y = k sin(29 B ) (0 B — 0)/\V(G)/ E\. y = gives 
the exact direction of the Bragg angle 6 = 8b- The crystal function vq describes the 
diffraction property of the wave function through a single crystal of thickness D. vq 
characterizes the case of diffraction due to the dynamical theory. The three terms in 
Eq.(l) form a superposition of three wave functions which belong to the three pathways 
(acf), (adg) and (beg). They have different phase shifts and absorptions. The crystal 
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functions defined above are equal for each of the three wave functions, because each beam 
is transmitted (function Vq) and diffracted (functions vq and vq(— y)) twice respectively. 

The wave function in the diffracted direction behind the analyzer crystal reads 

Here, ip e i = Mo e * kGr is a plane wave function with wave vector pointing to the 
diffracted direction. e lS is a phase factor. The crystal function v' G is equal to v (—y). 
From Eq. (4) it can be recognized that the three wave functions which are superimposed 
in the forward direction (0) have the same but those in the diffracted direction G do not 
have the same number of reflections and transmissions. Therefore an unsymmetrical 
intensity behavior can be expected in the diffracted direction, as shown below. 

In order to compute the intensities K and Kq behind the analyzer crystal the squared 
moduli of the wave functions, i.e. |^/)+(^)+(^) |2 and |^«c/)+(«&)+(fe ff ) ^ have to 

be taken into account. Because of the thickness of about 0.5 cm of the crystal plates 
in Fig.(l) the quantity A = nD/A » 1 cause very rapid oscillations of the terms 
sin 2n (x) which appear in these expressions. Therefore the mean values of these trigono- 
metric functions can be used. One gets: sin 2 (x) = 1/2, sin i (x) = 3/8, sin 6 (x) =5/16 
and sin s (x) = 35/128 [2]. Moreover integration over the variable y must be performed 
because of the beam divergence which has to be taken into account. Finally a normalized 
Gaussian spectral distribution 

g(k) = -Jr exp ^ k ~^f ] (5) 

of incoming wave numbers k is assumed. Here ko is the mean value and (5k) 2 is the 
mean square deviation of wave numbers. The intensities can now be calculated (setting 
Uq = 1) as follows (the limits of integration are always — oo and +oo): 

K = J 9(k)[ J \4 acf)+iad9)+ibe9) \ 2 dy}dk = ^{1 + e~ 2 ^ + e~ 2 ^ + 

+2e- a «e-W 2 ^ 2 / 2 cos(A d k ) + 2e~ a ' e"^) 2 ^/) 5 l / 2 cos(A f k ) + 
+2e -K+«/) e -W 2 (A d -A / )V2 cos[(Ad _ A f )k }} , (6) 

k g = J g(k)[j |4 ac/)+(a<i9)+(6e9) l 2 ^]^ = + ^ 2ad ) + 

+417e~ 2a ' + 130e- ad e" (<5A;)2(Ad)2/2 cos(A (i A;o) - 
-158e-^e-^ 2 ( A ^ 2 / 2 cos(A / A;o) - 

-158e' iad+a ^e- {Sk)2{Ad - A f )2/2 cos[(A d -A f )k }} . (7) 

Fig. (2) plots Kq and Kq for A^ = and a d = ctf = as a function of Af. If Af = 0, 
Kq = 71l7r/2048 and Kq = 36l7r/2048. These results can be found elsewhere [2]. One 
can recognize that the mean intensities of the two beams are - in general - different (just 
like in single loop interferometry) . The question arises if these intensity levels can be 
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Figure 2. Intensities Kq and Kq of Eqs.(6) and (7) behind the double loop 
interferometer (s. Fig.(l)) for a d = a f = 0, A d = and (6k) /ko = 1/100. The 
level of oscillations of Kq is below that of Kg- 



made equal by using special parameter values. The answer is yes. If one takes the mean 
values of K and K G to be equal, i.e. [K — K G ] Af = 0, one gets a condition for af 

01 f = ~\ ln{ Xm 1 1 + e ' 2ad + 2 e ~ ade ~ m2(Ad)2cos ( A dk )]}>o . (8) 

Any absorption a has to be >0, because the transmission T = e~ 2a of a beam obeys 
the relation 0<T<1. In case of a d = and A d = 0, the absorption af becomes 
— (l/2)/n(28/169) = 0.8988. Fig. (3) demonstrates this example. In this case the double 
loop interferometer may be regarded as a 50/50 beam splitter by using an absorption 
element in a suitable arm of the device. This is not possible using a single loop 
interferometer. 



3. Visibilities of single and double loop interferometers 

The visibility is defined as 

1 — TT~ • ( ) 

J-max ~ J-min 

We distinguish two kinds of absorption processes: stochastic (sto) and deterministic 
(det) absorption [6]. In an interferometer stochastic absorption is realized by inserting 
an absorbing material in one arm of a beam (s. Fig.(l)). Deterministic absorption 
can be achieved by a chopper which blocks the beam path periodically or by partial 
reduction of the beam cross section. In the following we investigate the visibility in a 
single and double loop interferometer at low interference order (coherence function is 
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Figure 3. Calculated intensity oscillations of K and Kq of Eqs.(6) and (7) 
considering the following parameters: a d — 0, aj = — (l/2)Zn(28/169) = O.i 
Ad = and (5k) /ko = 1/100. The intensity levels are equal due to a/. 



approximately 1). 

a) Single loop interferometer (index 1): If one phase shifter A d and one absorber a d 
is put in beam path (d), the intensity I in the forward direction can be written as 

I ^\l + e tk ° A ^\ 2 = l + T d + 2^T d cos(A d k ) , (10) 

where T d = e~ 2ad is the transmission probability of beam (d) and A d ko is the phase 
difference between the two beam paths. The visibility is therefore 

V stol = ^r ■ (11) 
J- + J-d 

For deterministic absorption we get 

J oc | 1 + e* ,Ad \ 2 T d + ( 1 - T d ) = 1 + T d + 2 T d cos(A d k ) , (12) 

W = T^f • (13) 
1 + J-d 

The difference between stochastic and deterministic absorption in a single loop neutron 
interferometer is drawn in Fig. (4) and has been discussed in [6]. 

b) Double loop interferometer (index 2): If A d and a d is placed in path (d) and no 
phase shifter or absorber in beam path (/) (see Fig.(l)), the intensity is given as (see 
Eq.(6)): 

K oc\2 + e tkaAd ~ ad \ 2 = T d + A[l + JT d cos(A d k )} , (14) 
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Figure 4. Visibilities of single and double loop interferometers. The visibilities V s to\ 
for stochastic (Eq.(ll)) and Vdeti for deterministic (Eq.(13)) absorption in a single 
loop interferometer are shown as a function of transmission probability T d in beam 
path (d) . In the double loop interferometer no phase shifter A f and no absorber a/ is 
in beam path (/). In this case the visibilities V st0 2 and Vdeti of Eqs.(15) and (17) are 
smaller than in a single loop device. 



V s to2 = < 1 • (15) 

4 + Id 

The deterministic case may be expressed as 

A'„(x|2 + e*> A ''| 2 r i + 4(l-r i ) =T d + 4[l + r d cos(A i *„)] , (16) 

V- = < 1 • (17) 

In general, V stol > V sto2 and V detl > V det2 for 0<T d <l (Fig.(4)). An interesting case 
arises, if a phase shifter Af (and no absorption «/) is inserted in beam path (/). The 
intensity then becomes 

K oc I e ik(>Ad - ad + 1 + e ikoA f\ 2 = 

= T d + 2 sJr d [cos{ A d k )+ cos( A d k -A f k )}+ 4 cos\ A f k Q /2) . (18) 

If Af = 0, then Eq.(18) reduces to Eq.(14). If A f k = (2n + 1)tt, then e iA f k ° = -1 
and K (xT d . If 2mr<Afko<(2n + l)ir, maxima of K are at A d k = Afk /2 + 2nn and 
minima at A d k = Afk /2 + (2n + l)ir. These maxima and minima are interchanged 
in the range (2n + l)vT<Aj/c <(2n + 2)n, and it is not necessary to consider this case 
separately. One gets: 

K ^<xT d ±4cos{A f k /2)[JT d ±cos(A f k /2)} , (19) 
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where the upper sign belongs to the maximum and the lower sign to the minimum 
intensity. The visibility is 

4VT d cos(A f k /2) 
Vsto2 ^- 4cos*(A f ko/2)+T d ■ (20) 

If Afk = 0, then V sto2 A f = V sto2 . If A f k = 2n/3, then V sto2 A f = V stol . In particular we 
emphasize the fact that the maximum of the function V st0 2A f is at T d = 4cos 2 (Aj/c /2) 
where V s t 2A f = 1- For a given transmission T d in beam path (d) a phase shift 
Afk = 2 arccos(\/Td/2) in path (/) should be chosen in order to achieve a visibil- 
ity of 1 (Fig. (5) and Fig. (6))! Setting T d — 1 the relation V st0 2A f of Eq.(20) reduces to 
a formula which already has been discussed in [3] (see Fig. (6)). 



An analogous equation can be derived for the deterministic case which reads (Fig. (7)) 

4T d cos{A f k /2) 
4cos 2 {A f k /2)+T d 

In case of deterministic visibility the value of 1 cannot be achieved for T d < 1. 



V det2A f = A _„„ A z /o^ i hp < 1 , for T d <l . (21) 



In order to compare intensities and to assess possibilities for measurements, Figs.(8a,b,c) 
present values of K and Kq as a function of the phase shift A d k$ in beam path (d). 
From these figures it can be concluded that in spite of strong absorption small signals 
can be measured with visibility 1 and with sufficient intensity by adequately adjusting 
the phase shifter in the second loop of the interferometer. However, as shown in the 
next section, the described procedure specifies an ideal situation and there is a limit for 
measuring weak signals using such a technique. 



4. Discussion and conclusion 



In this paper some intensity and visibility aspects of a double loop neutron interferom- 
eter have been considered. In the first part the wave functions of a triple wave function 
superposition in the forward and diffracted direction behind the interferometer has been 
described using the theory of dynamical diffraction in single crystals. Phase shifters and 
absorbers have been taken into account. In general the intensity levels of the forward 
and diffracted beams are different. Under certain circumstances these levels can be 
made equal using various absorbers. This is an interesting new aspect because in this 
case the device can be regarded as a 50/50 beam splitter. This feature is not possible 
in a single loop interferometer. 

In the second part the visibility has been investigated. As an another important new 
result it has been shown that a visibility value of 1 may be reached when the transmis- 
sion T d (for stochastic absorption processes) in the first loop is adjusted accordingly to 
the phase shift Afk in the second loop: 

T d = Acos 2 (A f k /2) -> V sto2Af = l. (22) 
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Figure 5. Visibility V s to2A f (Eq.(20)) in a double loop interferometer as a function of 
transmission probability Td in beam path (d). A phase shift Ayfco is applied in beam 
path (/). The visibility attains the value 1 for Td = 4 cos 2 (A/fc /2). 




Figure 6. Visibility V s to2A } (Eq.(20)) in a double loop interferometer as a function 
of phase shift Afk a in beam path (/). A transmission probability Td is applied in 
beam path (d). The visibility attains the value 1 for Td = 4cos 2 (A/fco/2), where the 
4 parameters T d correspond to those 4 parameters Afk in Fig. (5) via this equation. 
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Figure 7. Visibility Vd e t2A f (Eq.(21)) in a double loop interferometer as a function 
of transmission probability Td in beam path (d). A phase shift A/fco is applied in 
beam path (/). Contrary to V st0 2A f this function Vd e t2& } cannot attain a value of 1 
for T d < 1. 
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Figure 8. Intensities Ko and i^c behind the double loop interferometer (see Fig.(l)) 
as a function of the phase shift A d fc in beam path (d). a) No absorption in beam 
(d) (Td — 1) and no phase shift in beam (/) (A/fco = 0). Fig. (8a) can be compared 
to Fig. (2). Note that in Fig. (8a) the intensities K and Kq are plotted as a function 
of Adko. b) Absorption in beam (d) (transmission Td = 0.1) and no phase shift in 
beam (/). The intensities have been attenuated accordingly, c) A/fco = 2.824 has 
been chosen. This value corresponds to Td = 0.1 (see Figs. (5) and (6)). K n shows a 
visibility of 1 because of I m in = 0. 
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Td = e~ 2ad is the transmission probability in beam (d) (a^ is the absorption coefficient) 
and Af is the phase shift in beam (/), where - according to Fig.(l) - the intensity K in 
the forward direction behind the double loop interferometer oscillates as a function of 
the phase shift A^ in beam (d). Using the described method in principle every signal in 
path (d) can be transformed to a signal Kq with visibility 1 and therefore the method 
could be most interesting for small signals. 

It should be noted that to describe real experimental situations an additive term Ii nco h 
must be introduced in Eq.(18) to represent the incoherent part of the intensity (back- 
ground intensity). The main reason for the incoherent effects is non - interference 
because of crystal imperfection (not absolute parallelism of the crystal lattice planes 
throughout the interferometer) as well as lattice vibrations and small temperature gra- 
dients. The term Ii nco h leads to a new visibility V' sto2Aj : 



This expression is always less than 1 for Ii nco h > (see also Eq.(22)). In order to mea- 
sure weak signals (see Fig. (8c)) the background intensity has to be as small as possible. 
This could be a serious constraint considering real experimental conditions using neu- 
trons. It should be mentioned that the background intensity has no influence on Eq.(8) 
concerning the 50/50 beam splitter system. Note also that an additional empty phase 
which is a signature of an interferometer does not affect the aforementioned considera- 
tions about the intensity and visibility. 

It is a common problem in signal processing to measure weak signals. The absorber 
serves only to generate these signals. In spite of the constraints mentioned above for 
almost ideal experimental situations Eq.(22) predicts a visibility close to 1. This could 
also be the case in systems of Mach - Zehnder interferometers using laser beams. Here 
the number of counts registered in a detector in a certain time interval obeys a Poisson 
distribution with variance (AN) 2 = N, where N is the mean number of counts. The 
Poisson distribution is a signature of coherent state behavior. For small count rates the 
described method may be used to measure small signals with visibility 1. 

The system can easily be adapted (loop B in Fig.(l)) to an eight-port interferome- 
ter as a high sensitive homodyne detection setup [4,7,8] as shown in Fig. (9). In this 
case two vacuum inputs |0 > have to be added and the strong local oscillator limit 
can be applied to the idler beam (a) when the signal beam (e) is relatively strongly 
attenuated by an absorber a e . Homodyne detection in laser interferometry is a method, 
in which the field amplitudes (the quadrature components) are measured instead of the 
quantized intensity. The balanced version of homodyne detection has a great practical 
advantage of cancelling technical noise and the classical instabilities of the reference 
field. Thereby the signal interferes with a coherent laser beam at a well-balanced 50/50 
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Figure 9. An eight-port interferometer is built up from loop B in Fig.(l). Neutron 
count rate differences can be measured by the detectors and Di as well as D fi and 
D g i . The phase of the strongly attenuated signal beam (e) can be modified by a phase 
shift A e . The strong local oscillator is given by the idler beam (a). 

beam splitter. It provides the phase reference for the quadrature measurement. The 
intensity difference is the quantity of interest because it contains the interference term 
of the local oscillator and the signal. 

Finally, as already mentioned above, we would like to point out again that there are 
good reasons to assume that a similar result can be attained in double Mach-Zehnder 
interferometry where photonic beams in fiber glass are used. As a visibility near to 1 
is of great advantage for measuring largely attenuated beams, our approach outlines a 
general method to investigate weak signals in double loop interferometric devices. 
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